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The  construction  of  an  optimum  model-sampling  simulation  experimental 
design  is  often  an  art  and  clearly  very  much  problem  dependent.  The  accurate 
estimation  of  the  distribution  of  any  random  variable  by  computer  is  almost 
always  difficult,  especially  when  reasonable  accuracy  is  desired  for  its  tail 
probabilities.  If,  in  addition,  each  element  of  the  computer  generated  pseudo¬ 
random  sequence  is  Itself  the  result  of  a  stochastic  limit  or  possibly  a 
functional  of  a  continuous-time  process,  it  becomes  most  complicated  to  assess 
the  final  statistics.  In  this  paper,  we  focus  on  the  estimation  of  tail 
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probabilities  for  the  distribution  function  of  the  maximum  on  the  unit 
interval  of  a  continuous-time  Wiener  process  approximated  as  a  multi¬ 
variate  normal  of  increasing  dimension.  We  critique  recent  approaches 
to  sample-size  determination  for  such  distribution-sampling  problems 
and  build  on  insights  from  probability  theory  to  find  more  reasonable 
run  sizes. 
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Abstract 


The  construction  of  an  optimum  model -sampl i ng  simul  ^n 
experimental  design  is  often  an  art  and  clearly  very  much  prc  iem 
dependent.  The  accurate  estimation  of  the  distribution  of  any 
random  variable  by  computer  is  almost  always  difficult, 
especially  when  reasonable  accuracy  is  desired  for  its  tail 
probabi 1 i ti es.  If,  in  addition,  each  element  of  the  computer, 
generated  pseudo-random  sequence  is  itself  the  result  of  a 
stochastic  limit  or  possibly  a  functional  of  a  continuous-time 
process,  it  becomes  most  complicated  to  assess  the  final 
statistics.  In  this  paper,  we  focus  on  the  estimation  of  tail 
probabi 1 i ti es  for  the  distribution  function  of  the  maximum  on 
the  unit  interval  of  a  continuous-time  Wiener  process 
approximated  as  a  multivariate  normal  of  increasing  dimension. 
We  critique  recent  approaches  to  sample-size  determination  for 
such  distribution-sampl ing  problems  and  build  on  insights  from 
probability  theory  to  find  more  reasonable  run  sizes. 


c _ 


1 •  Introduction 

The  determination  of  a  best  run  size  for  a  model -sampl i ng 
simulator  is  generally  a  nontrivial  matter.  Many  authors  (for 
example,  see  Bratley,  Fox  and  Schrage,  1983,  and  Rubinstein, 
1981)  have  commented  on  the  difficulty  of  accurately  estimating 
measures  of  random  variables  using  computers  even  when  generating 
random  (or  pseudo-random)  samples.  Limit  theorems  are  typically 


of  little  use  when  the  underlying  probabilistic  structure  is 
complex  and  rates  of  convergence  not  available.  Furthermore,  it 
is  often  quite  hard  to  obtain  appropriate  estimates  of 
variability  for  calculating  confidence  intervals  or  making 
statements  of  accuracy. 

It  is  well  known  that  run  lengths  necessary  to  achieve 
desired  levels  of  "precision"  for  estimating  such  things  as 
population  moments  and  quantiles  may  rise  very  quickly  with  the 
order  of  the  percentile  or  moment.  With  the  advent  of 
inexpensive  and  accessible  computing,  resource  and  time 
constraints  are  generally  not  at  issue  (though  repeated  sorting 
may  be  necessary).  Instead,  it  often  becomes  more  important  to 
understand  the  behavior  of  the  population  parameters  sought. 
This  sort  of  problem  becomes  even  more  complicated,  when  the 
values  of  each  element  of  an  I ID  sequence  are  themselves  the 
result  of  a  further  limiting  process,  (as,  for  example,  in 
extreme-value  problems),  so  that  the  full  sequence  is  actually 
one  small  series  imbedded  within  another  larger  one. 

Excellent  examples  of  these  problems  arise  in  the  use  of 
simulation  methods  for  estimating  functionals  of  Gaussian 
processes.  We  draw  attention  to  the  work  of  Serfling  and  Wood 
(1976),  Wood  (1978),  Siegmund  (1978),  and  Chandra,  Singpurwalla 
and  Stephens  (1983).  The  consi derati on  of  such  problems  is 
important  in  statistics  primarily  because  of  their  appearance  in 
the  asymptotic  theory  of  goodness-of-f i t  tests  (see  Chandra, 
Singpurwalla  and  Stephens,  1980  and  1981). 

One  of  the  highlighted  illustrations  from  the  1983  Chandra, 
Singpurwalla  and  Stephens  paper  (henceforth  called  CSS)  offers  an 


excellent  example  of  the  difficult  search  for  a  satisfactory 
experimental  design  and  ultimately  the  need  to  make  very  long 


computer  runs  for  adequate  precision.  This  problem  is  the 
estimation  by  simulation  of  upper  quantiles  for  the  maximum  of 
the  (continuous-time)  Brownian  motion  process  over  the  unit 
interval  (call  it  W(t),  0  <  t  <  1).  We  know  that 


ECW(t)]  =  0  (for  all  t;  W(0>  =  O) 

1 

and  that 

EEW(t)W(s)3  =  min(s,t)  (0  <  s,t  <_  1).  (1) 

^  From  a  practice  first  used  by  Serf  ling  and  Wood,  we  can 

approximate  this  process  at  the  points  (j/k,  k  fixed  and  O  j 
<_  k)  by  (k+1 ) -vari ate  normal  vectors  with  mean  0  and  the  same 
covariances  as  the  Brownian  motion.  We  know  that 

I 

lim  Prfmax  W(j/k)  >  xJ  =  2<l-®(x)>  (x  >0),  (2) 

k-s**  j<k 

where  ®(x)  is  the  univariate  normal  CDF.  Furthermore,  Siegmund 
showed  that  the  complementary  CDF  of  the  maximum  could  be  well 
approximated  for  finite  k  as 

! 

i 

I 

t 

;  2C1-® (x+. 583/ Vk) >  ,  (3) 


with  the  result  binding  in  the  limit  as  k  ■*  *.  Note  also  that 


the  limiting  density  is  the  -folded  standard  unit  normal  on 

We  often  see  estimation  of  both  the  population  mean  and 
variance  from  the  same  pseudo-random  sample.  However,  the 
variance  of  a  percentile  is  itself  a  function  of  the  probability 
distribution  of  the  variable.  More  completely,  we  know  (for 
example,  see  Durbin,  1973)  that  the  properly  scaled  and 
translated  (np)th  order  statistic  from  an  ordinary  random  sample 
has  a  limiting  normal  distribution,  that  is,  that 

*  (np)  ""^p  7 

-  — >  Z  ~  N <0 , 1 /f " < £p ) )  (f  =  original  density). 

Vp ( 1-p) /n 

Since  this  variance  of  the  upper  quantiles  depends  on  tail 
probabi 1 ities,  it  is  clear  that  run  size  must  rise  rapidly  as  p 
goes  to  one  for  decaying  densities.  For  example,  the  unit  normal 
has  approximate  functional  values  of  .1755  and  .02665  for  p=.9 
and  .99,  respect i vel y .  Hence  the  estimator  of  the  (,99n)  order 
statistic  has  an  approximate  variance  of  13.94/n  as  compared  to 
2.922/n  for  the  (,9n)  one,  or  a  nearly  five-fold  increase.  This 
further  reinforces  our  earlier  point  of  the  need  to  make  longer 
runs  for  the  fixed  precision  estimation  of  tail  quantiles. 

Our  target  in  this  work  then  was  to  develop  a  more  complete 
perspective  on  appropriate  run  sizes  for  this  Gaussian  problem. 
Of  course,  the  exact  critical  values  are  known  here;  but  this 
affords  an  excellent  opportunity  to  calibrate  any  approach  for 
determining  sample  sizes. 


Resul ts 
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CSS  generated  10,000  ( k+1 ) -var i ate  normal  vectors 
(k=20,30,50,60,90)  using  the  extended  precision  version  of  the 
routine  GGNSM  -from  the  International  Mathematical  and  Statistical 
Library  (IMSL).  Their  Table  I  shows  the  pttn  quantiles  of  the 
empirical  distribution  of 

max  fW ( j/k) > 

j<k 

from  their  simulation  experiment  for  p  =.9,  .95,  .975  and  .99, 
using  10,000  replications.  These  results  are  compared  to  the 
correspondi ng  quantiles  given  by  Siegmund's  approx i mat i on ,  and 
the  authors  claim  to  have  encountered  an  important  anomaly  in 
such  a  simulation,  namely,  that  there  is  a  critical  "turning  down 
effect  as  k  becomes  large,"  which  puts  their  point  estimates  on  a 
path  below  what  they  should  be.  Even  if  we  assume  that  the 
Siegmund  approximation  is  close  to  correct  for  smaller  k,  there 
should  be  some  variability  about  its  values  for  any  run  size  by 
the  very  nature  of  random  draws.  However,  CSS  felt  that  they  had 
seen  too  much  movement,  and  even  a  pattern  down  and  away.  But  we 
observe  that  the  kind  of  variation  they  noted  is  totally  to  be 
expected  in  light  of  what  is  really  a  rather  small  sample  size, 
and  the  possible  need  for  vastly  different  run  sizes  for 
different  p  values. 

In  our  mind  the  problem  of  determining  rates  of  convergence 
of  sample  quantiles  for  the  maximum  of  a  Brownian  motion  by 
simulation  is  very  challenging  and  a  good  illustration  for  a  more 


•V  *.  V 


general  class  of  such  problems.  So  we  attempted  to  examine  the 
problem  more  carefully  and  to  develop  a  sounder  strategy  for 
finding  an  agreeable  run  size.  We  did  this  by  repeating  the 
basic  CSS  experimental  design.  We  have  created  a  F0RTRAN77 
program  using  GGNSM,  with  the  same  values  of  k ,  namely, 
20,  30,  50,  60  and  90,  and  with  p-values  of  .9,  .95,  .97  and  .99. 
(We  have  opted  to  use  .97  instead  of  .975  to  keep  the  abscissa 
spacing  uniform).  Complete  runs  were  performed  for  the  eight 
sample  sizes  of  10,000,  20,000,  30,000,  40,000,  50,000,  70,000, 

90 , 000  and  1 00 , 000 . 

Our  simulation  results  are  all  documented  in  Tables  I -VI II, 
where  they  are  compared  to  those  of  CSS  and  what  might  be 
expected  from  the  Siegmund  approximation.  To  put  all  these 
numbers  in  clearer  perspective,  we  have  made  two  sets  of  plots. 
The  first  set  is  represented  by  Figures  1-3,  where  the  quantiles 
for  the  three  p  values  are  plotted  for  the  largest  k  =  90  and 

compared  to  their  approx i mat i on  given  by  (3).  The  second  group 
provides  a  picture  of  the  estimated  rate  of  convergence  of  the 
critical  values  as  k  increases.  As  done  by  CSS,  on  each  of  these 
figures  we  have  included  a  plot  of  the  line  which  results  for 
each  quantile  from  Equation  (3). 

To  summarize  the  salient  insights  offered  by  the  first  three 
figures,  we  note  that  the  results  are  not  at  all  surprising.  For 
the  moment  assuming  that  the  exact  vari abi 1 i ti es  of  the 
estimators  are  not  known,  each  estimate  (i.e.,  the  .90,  .95  and 
.99  quantiles)  has  become  moderately  stable  with  the  increasing 
sample  sizes  (for  k  =  90).  As  expected,  the  .99  estimate  is  the 


least  stable  of  the  three.  We  might  even  find  all  of  the  100,000 


sample  values  acceptably  precise  to  two  or  even  three  significant 
digits  (but  certainly  no  more  than  that). 


However,  this  conclusion  may  be  very  misleading.  A  more 
completely  constructed  experimental  design,  arranged  into  blocks 
for  estimating  variances,  would  likely  show  a  rather  wide 
confidence  interval  for  the  .99  quantile,  down  to  a  fairly  tight 
one  for  the  .90th.  In  light  of  the  fact  that  the  underlying 
distribution  theory  is  available  for  the  Brownian  maximum  (that 
is,  from  Equation  (2)),  we  can  actually  compute  the  variance  of 
each  of  our  three  key  quantiles.  Since  the  density  function  of 
the  maximum  is  the  folded  unit  normal,  we  are  easily  able  to 
compute  the  appropriate  tail  ordinates  as 


.  350996 

for 

p  =  .90 

. 206272 

for 

p  =  .95 

m 

. 053304 

for 

p  =  .99 

Thus  it  follows  that  the  standard  deviations  for  the  estimates  of 
the  .90,  .95  and  .99  quantiles  taken  from  a  sample  of  100,000  are 
found  as 


p (1-p) 

nf2(£p> 


!. 0027028 
. 0033412 
. 0055903 


(p  =  .90) 
(p  =  .95) 
(p  =  .99) 


For  purposes  of  understandi ng  the  precision  of  our 
estimators,  let  us  overlay  these  (limiting)  standard  deviations 
onto  the  k  =  90  final  estimates  of  Figures  1-3  to  make 
approximate  (say,  95%)  confidence  statements.  Thus  the  following 
intervals  result  (compared  to  Siegmund's  estimate  for  k  =  90): 


(1.5893,1.6001)  -for  p  =  .90  vs.  Si  egmund  '  s  value  of  1.5835 
(1.9026,1.9160)  for  p  =  . 95  vs.  Si egmund 's  value  of  1.8985 
(2.5174,2. 5398 )  for  p  =  .99  vs.  Siegmund's  value  of  2.5385 


Of  course,  we  note  the  much  wider  interval  for  .99.  In  fact,  we 
can  calculate  the  sample  size  which  would  be  needed  to  give  the 
same  absolute  accuracy  for  the  .99  estimate  as  in  .90.  This 
would  be 


N 


(.99)  (.01) 

- - - -  *=  477,000, 

( . 0027028) 2 ( . 053304 ) 2 


or  nearly  f i ve  times  the  baseline  sample. 

For  a  sample  of  10,000,  the  standard  deviations  for  the 
quantile  estimates  increase  V10-fold  to 


p ( 1 -p ) 
nf2(£p) 


{.0085471  ( p  = . 90 ) 

.0105659  (p= . 95 ) 

.0176781  (p=. 99) 


The  resultant  confidence  intervals  (for  any  level)  must  thus  be 
VI 0  times  as  wide.  Thus,  for  example  based  on  k=90  the  actual 
.99  point  would  be  estimated  to  be  within  (2. 5428-. 0353562, 
2. 5428+. 0353562)  or  (2.5074,  2.5782),  which  is  quite  a  broad 

interval.  In  our  view,  therefore,  10,000  is  just  too  small  a  run 
size  to  give  adequate  precision,  or  indeed  maybe  to  prevent  any 
conclusion  whatsoever. 

The  primary  lessons  of  Figures  4-6  focus  on  the  rates  of 


convergence  with  respect  to  the  parameter  k 


The  work  of 


Siegmund  clearly  suggested  that  larger  values  of  k  should  be  used 


■for  more  accurate  estimation.  But  there  is  obviously  a  -fairly 
large  increase  in  computing  time  as  the  dimension  o-f  the 
multivariate  normal  goes  up;  but  the  payo-f-f  in  precision  could  be 
significant,  since  simple  linear  extrapolation  may  not  work  well. 

The  result  o-f  Equation  (3)  translates  into  a  straight  line 
-for  each  quantile  when  plotted  against  1/Vk=x.  For  .90  the  line 
is  y=  1 . 645—.  583x  ,  while  -for  .99  it  is  y=2.  598-.  583>: .  Since  the 
slopes  are  the  same  for  all  quantiles,  each  line  reaches  the  y— 
axis  at  a  point  .583/V90  -  .0614  higher  than  that  for  x  =  l/-/90. 
Thus,  for  example,  to  reach  a  point  within  .01  of  the  limiting 
height  would  require  a  k  of  58.  3^  or  about  3,400.  Furthermore, 
we  notice  that  all  95"/.  confidence  intervals  constructed  for  k=90 
do  not  cover  the  actual  answers.  Hence,  it  seems  clear  that 
extreme  caution  should  be  used  with  this  sort  of  approach. 

To  combine  the  results  displayed  in  the  figures,  we  might 
conclude  that  the  .90  quantile  is  reasonably  well  approximated 
with  a  sample  size  of  100,000  and  k  =  90,  though  a  somewhat 
larger  k  would  be  even  better.  For  the  .95  quantile,  a  slightly 
larger  run  size  would  be  preferable,  with  k  definitely  increased 
beyond  90.  And,  finally,  to  get  adequate  precision  for  the  .99 
point,  even  larger  increases  would  be  warranted  in  the  run  size 
and  k  possibly  up  to  as  much  as  a  sample  of  500,000  and  an  order 
of  magnitude  (or  higher)  increase  in  k. 
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TABLE  I 


Estimated  Percentage  Points  for  Run  Size  =  10,000 


Values 

Values 

CSS 

From 

Harris 

of  k 

o-f  p  Estimate 

Approx 

Estimate 

20 

.90 

1 . 5725 

1.5146 

1.5367 

.95 

1.9318 

1 . 8296 

1.8193 

.97 

- 

# 

2.  0275 

.99 

2.5053 

2.4696 

2.  4797 

30 

.90 

1 . 6049 

1 . 5386 

1 , 5560 

.95 

1 . 9362 

1 . 8536 

1.8571 

.97 

- 

* 

2.  0592 

.99 

2.5815 

2.4936 

2.4536 

50 

.90 

1 . 6376 

1 . 5626 

1 . 5938 

.95 

1 . 9798 

1 . 8776 

1 . 8994 

.97 

- 

* 

2.0697 

.99 

2.5407 

2.5176 

2.4565 

60 

.90 

1 . 5997 

1 . 5697 

1.5712 

.95 

1.9416 

1 . 8847 

1 . 8767 

.97 

- 

* 

2.0865 

.99 

2.4642 

2.5250 

2.5264 

90 

.90 

1 . 6239 

1 . 5835 

1 . 5893 

.95 

1.8913 

1 . 8985 

1 . 8985 

.97 

- 

* 

2. 1189 

.99 

2.4625 

2.5385 

2.5428 

NOTE: 

The  hyphens 

are  in 

the 

.  97  rows  because  CSS 

did  not  estimate  -for 

that 

p  value,  while 

we  felt 

that  such  an 

estimate  provided  important 

insight. 

The  asterisk! 

s  are 

used 

to  indicate 

that  the 

approx i mat i on 

can  be 

used 

for  .97,  but 

we  chose 

not  to  do  the 

calculation 

since  a  full  comparison 

is  not  passible. 


TABLE  II 


Estimated  Percentage  Points  -for 

Run  Size  = 

20 , 000 

Values 

Values  CSS 

From 

Harri s 

of  k 

of  p  Estimate 

Appro* 

Estimate 

20 

.90 

1 . 5725 

1.5146 

1 . 5352 

.95 

1.9318 

1 . 8296 

1 . 8429 

.97 

- 

* 

2.0525 

.99 

2. 5053 

2.4696 

2.4700 

30 

.90 

1 . 6049 

1 . 5386 

1.5271 

.95 

1 . 9362 

1 . 8536 

1 . 8402 

.97 

- 

* 

2.0760 

.99 

2.5815 

2. 4936 

2.4918 

50 

.90 

1 . 6376 

1 . 5626 

1 . 5594 

.95 

1 . 9798 

1 . 8776 

1.8612 

.97 

- 

* 

2. 1042 

.99 

2.5407 

2.5176 

2.5044 

60 

.90 

1 . 5997 

1 . 5697 

1 . 5533 

.95 

1.9416 

1 . 8847 

1 . 8599 

.97 

- 

* 

2.0954 

.99 

2.4642 

2.5250 

2.5102 

90 

.90 

1.6239 

1 . 5835 

1 . 5838 

.95 

1.8913 

1 . 8985 

1 . 8777 

.97 

- 

2. 1885 

2.0691 

.99 

2.4625 

2.5385 

2.4787 

TABLE  III 


Estimated  Percentage  Points  -for  Run  Size  =  30 


Values  CSS  From 

o-f  p  Estimate  Approx 


90 

1 . 5725 

1.5146 

95 

1.9318 

1 . 8296 

97 

- 

* 

99 

2.5053 

2.4696 

90 

1 . 6049 

1 . 5386 

95 

1 . 9362 

1 . 8536 

97 

- 

# 

99 

2.5815 

2. 4936 

90 

1 . 6376 

1 . 5626 

95 

1 . 9798 

1 . 8776 

97 

- 

* 

99 

2.5407 

2.5176 

90 

1 . 5997 

1 . 5697 

95 

1.9416 

1 . 8847 

97 

- 

# 

99 

2.4642 

2.5250 

.90 

.95 

.97 

.99 


1.6239 

1.8913 


1 . 5835 
1 . 8985 
2. 1885 
2.5385 


000 


Harris 

Estimate 


1.5312 

1.8419 

2.0541 

2.4679 

1.5361 
1 . 8584 
2.0505 
2. 4630 

1.5613 

1.8731 

2.0742 

2.4971 

1 . 5690 
1 . 8923 
2. 0932 
2.4837 

1 . 5666 
1 . 8892 
2. 0927 
2.5068 


2.4625 


TABLE  IV 


Values 
o-f  k 

20 

30 

50 

60 

90 


Estimated  Percentage  Points  -for  Run  Size  =  40 


Values  CSS  From 

o-f  p  Estimate  Approx 


90 

1.5725 

1.5146 

95 

1.9318 

1 . 8296 

97 

- 

* 

99 

2.5053 

2.4696 

90 

1 . 6049 

1 . 5386 

95 

1 . 9362 

1 . 8536 

97 

- 

* 

99 

2.5815 

2. 4936 

90 

1 . 6376 

1 . 5626 

95 

1 . 9798 

1 . 8776 

97 

- 

* 

99 

2.5407 

2.5176 

90 

1.5997 

1 . 5697 

95 

1.9416 

1 . 8847 

97 

- 

* 

99 

2.4642 

2.5250 

.90 
.95 
.  97 
.99 


1.6239 

1.8913 


1 . 5835 
1 . 8985 
2. 1885 
2.5385 


000 


Harri s 
Estimate 


1.5281 

1 . 8356 
2. 0505 
2.4679 

1 . 5375 
1 . 8404 
2.0514 
2.5039 

1 . 5594 
1 . 8707 
2.0739 
2.4918 

1 . 5666 
1 . 8768 
2. 0837 
2.5144 

1 . 5846 
1 . 8977 
2. 1053 
2.5082 


2.4625 


TABLE  V 


Estimated  Percentage  Points  for  Run  Size  =  50,000 


Values 
of  k 


Values  CSS  From 

of  p  Estimate  Approx 


Harris 

Estimate 


20 

.90 

1 . 5725 

1.5146 

1 . 5284 

.95 

1.9318 

1 . 8296 

1 . 8358 

.97 

— 

* 

2.0505 

.99 

2.5053 

2.4696 

2.4648 

30 

.90 

1 . 6049 

1 . 5386 

1 . 5304 

.95 

1 . 9362 

1 . 8536 

1.8441 

.97 

— 

* 

2.0660 

.99 

2.5815 

2. 4936 

2. 4595 

50 

.90 

1 . 6376 

1 . 5626 

1 . 5667 

.95 

1 . 9798 

1 . 8776 

1.8781 

.97 

- 

* 

2.0801 

.99 

2.5407 

2.5176 

2.4777 

60 

.90 

1 . 5997 

1 . 5697 

1 . 5638 

.95 

1.9416 

1 . 8847 

1 . 8832 

.97 

— 

* 

2. 0982 

.99 

2.4642 

2.5250 

2. 4852 

90 


.90 

.95 

.97 

.99 


1 . 6239 
1.8913 


1 . 5835 
1 . 8985 
* 

2.5385 


1 . 5797 
1 . 8894 
2.0933 
2.4919 


2.4625 


TABLE  VI 


Estimated  Percentage  Points  for 

Run  Size  = 

70,000 

Values 
of  k 

Values  CSS 

of  p  Estimate 

From 

Approx 

Harri s 
Estimate 

20 

.90 

1 . 5725 

1.5146 

1 . 5264 

.95 

1.9318 

1 . 8296 

1 . 8338 

.97 

- 

* 

2.0505 

.99 

2. 5053 

2.4696 

2.4547 

30 

.90 

1 . 6049 

1 . 5386 

1.5410 

.95 

1 . 9362 

1 . 8536 

1.8501 

.97 

- 

* 

2. 0608 

.99 

2.5815 

2. 4936 

2.4731 

50 

.90 

1.6376 

1 . 5626 

1 . 5620 

.95 

1 . 9798 

1 . B776 

1 . 8857 

.97 

- 

* 

2.0895 

.99 

2. 5407 

2.5176 

2.4772 

60 

.90 

1 . 5997 

1 . 5697 

1 . 5704 

.95 

1.9416 

1 . 8847 

1 . 8878 

.97 

- 

* 

2.0921 

.99 

2.4642 

2.5250 

2.4996 

90 

.90 

1 . 6239 

1 . 5835 

1 . 5972 

.95 

1.8913 

1 . 8985 

1 . 9099 

.97 

- 

* 

2. 1177 

.99 

2.4625 

2. 5385 

2.5283 

-15 


TABLE  VII 


Val ues 
of  k 

20 

30 

50 

60 

90 


Estimated  Percentage  Points  for  Run  Size  =  90 


Values  CSS  From 

of  p  Estimate  Approx 


90 

1 . 5725 

1.5146 

95 

1.9318 

1 . 8296 

97 

- 

# 

99 

2.5053 

2.4696 

90 

1 . 6049 

1 . 5386 

95 

1 . 9362 

1 . 8536 

97 

- 

* 

99 

2.5815 

2.4936 

90 

1 . 6376 

1 . 5626 

95 

1 . 9798 

1 . 8776 

97 

- 

* 

99 

2.5407 

2.5176 

90 

1 . 5997 

1 . 5697 

95 

1.9416 

1 . 8847 

97 

- 

* 

99 

2.4642 

2. 5250 

.90 

.95 

.97 

.99 


1 . 6239 
1.S913 


1 . 5835 
1 . 8985 
# 

2.5385 


000 


Harri s 
Estimate 


1 . 5257 
1 . 8324 
2. 0505 
2. 4582 

1.5417 
1 . 8562 
2.0656 
2. 4732 

1 . 5688 
1 . 8856 
2. 0966 
2.4940 

1 . 5695 
1 . 8909 
2. 1056 
2.5153 

1 . 5960 
1.9127 
2. 1312 
2.5344 


2.4625 


TABLE  VIII 


Estimated  Percentage  Points  for  Run  Size  =  100,000 


Values 

Values 

CSS 

From 

Harr i s 

of  k 

of  p 

Estimate 

Approx 

Estimate 

20 

.90 

1 . 5725 

1.5146 

1 . 5240 

.95 

1.9318 

1 . 8296 

1 . 8320 

.97 

- 

* 

2.0492 

.99 

2.5053 

2.4696 

2.4565 

30 

.90 

1.6049 

1 . 5386 

1 . 5437 

.95 

1 . 9362 

1 . 8536 

1.8613 

.97 

- 

* 

2.0695 

.99 

2.5815 

2.4936 

2.4691 

50 

.90 

1 . 6376 

1 . 5626 

1 . 5646 

.95 

1 . 9798 

1 . B776 

1 . 8796 

.97 

— 

* 

2.0906 

.99 

2. 5407 

2.5176 

2. 4873 

60 

.90 

1 . 5997 

1 . 5697 

1 . 5774 

.95 

1.9416 

1 . 8847 

1 . 9037 

.97 

- 

* 

2.1120 

.99 

2.4642 

2.5250 

2.5192 

90 

.90 

1 . 6239 

1 . 5835 

1 . 5947 

.95 

1.8913 

1 . 8985 

1 . 9093 

.97 

— 

* 

2.1193 

.99 

2.4625 

2.5385 

2.5286 

99  QUANTILE  FOR  MAX  OF  BROWNIAN  MOTION  (K=90) 


20.00  30.00  40.00  50.00  60.00  7 0.00  80.00  90.00  100.00 

SAMPLE  SIZE  (IN  THOUSANDS) 
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